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Summary 

Conductors  carrying  very  high  currents  show  losses 
of  an  electromagnetic  and  a  shock  compression  nature. 
Electromagnetic  losses  (joule  heating  in  the  skin  lav- 

O  O  ° 

er,  magnetic  flux  diffusion)  scale  as  H  or  (I/a)  , 
where  H  is  the  self  magnetic  field  of  the  current  and 
I/a  is  the  current  1- conductor  periphery;  shock  losses 

scale  as  or  |(l/a)**  or  (l/a)^f  depending  on  the 
magnitude  of  I.  In  experiments  where  electrical  energy 
must  be  converged  from  a  large  pulsed  power  supply  to  a 
small  load,  these  losses  can  account  for  half  the  orig- 
nal  energy  and  limit  the  magnitude  of  the  energy  per 
unit  volume  in  the  load.  These  considerations  may  be 
important  for  the  study  of  material  properties  at  high 
energy  density. 

The  properties  of  materials  in  a  state  of  high  en¬ 
ergy  density  if  of  interest  to  a  broad  range  of  scien¬ 
tific  disciplines  ranging  from  astrophysics  to  material 
science.  Recent  advances  in  pulsed  power  technology 
have  made  possible  systems  with  energy  storage  in  the 
megajoule  range  and  have  made  possible  the  production 

of  material  energy  densities  ~  1  MJ/cm  . 

A  traditional  approach  to  achieving  high  energy 
density  is  to  obtain  such  a  high  energy  pulsed  power 
supply  and  design  a  "load"  of  arbitrarily  small  volume 
in  order  to  attain  a  desired  energy  density.  Where  the 
energy  store  is  in  the  form  of  charge  or  inductive 
storage,  the  load  is  connected  to  the  energy  storage 
device  by  a  suitable  transmission  line;  and  the  size  of 
the  line  and  the  load  are  often  constrained  by  mechan¬ 
ical  and  practical  electrical  considerations. 

Since  megajoule  class  energy  storage  systems  are 
very  large  —  many  meters  in  dimension  —  and  the  load 
is  often  measured  in  centimeters,  it  is  necessary  to 
devise  transmission  lines  which  converge  the  flow  of 
current  from  the  source  so  that  the  electrical  energy 
can  be  fed  into  the  load.  Near  and  at  the  load,  the 
power  convergence  is  maximum.  It  is  the  purpose  of 
this  paper  to  show  that  under  such  conditions,  electro¬ 
magnetic  and  shock  compression  losses  in  the  conductors 
can  account  for  a  significant  fraction  of  the  total  en¬ 
ergy  available.  The  losses  scale  as  (i/a)^  and  (i/a)1*, 
where  I/a  is  the  current  per  unit  length  of  current 
front  (e.g.  ,  for  a  circular  conductor,  l/a  =  current  4- 
circumference) .  Since  I/a  is  a  function  of  the  load 
geometry,  there  may  exist  an  optimum  lead  size  for  a 


maximum  energy  density  in  the  load.  Thus,  the  highest 
energy  density  is  not  produced  by  the  smallest  load 
volume.  The  presence  of  a  current-carrying  plasma 
sheet  on  the  cathode  side  of  the  transmission  line  may 
produce  a  modest  reduction  in  the  magnitude  of  the 
loss,  but  should  not  affect  the  scaling  with  (l/a). 

The  losses  we  consider  here  derive  from  the  pene¬ 
tration  of  the  self  magnetic  field  into  the  conductors 
as  well  as  the  associated  ohmic  heating  of  the  current- 
carrying  layer  in  the  conductors.  These  we  identify  as 
the  electromagnetic  energy  losses.  We  also  consider 
losses  due  to  the  shock  compression  of  the  conductor 
surface  layers  (the  shock  is  driven  by  the  external 
magnetic  pressure).  Unless  otherwise  stated,  the  units 
are  practical  CGS.  Unless  otherwise  stated,  references 
to  Knoepfel  are  to  "Pulsed  High  Magnetic  Fields,"  by 
H.  Knoepfel  (North  Holland  Publishing  Company,  1970). 


A.  Electro  Magnetic  Energy  Losses 

From  Poynting's  theorem,  the  total  electromagnetic 
energy  per  unit  area  flowing  into  the  conductor  (planar 
geometry)  is 
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=  -J±  J- 
4  (4T  T 


fgl  H(oji)dt 

J  2*  k-t 


VV-  10 

47T 

To  obtain  8H  (o,t),  one  needs  to  solve  the  magnetic 
dX 

diffusion  equation, 

•  kb~- 


The  solution  depends  on  the  time  dependance  of  the  mag¬ 
netic  field  at  the  boundary  surface.  Generally,  the 
pulse  shape  in  inductance-dominated  systems  approxi¬ 
mates  a  half-sinusoid;  we  approximate  this  with  a  time 
dependence  given  by 


The  integration  over  time  implied  in  the  energy 

equations  would  be  over  the  interval  o  <  t  <  t  .  For 

m 

the  proposed  form  of  H  (o,t),  Knoepfel  gives  the  so¬ 
lution  as 
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This  analysis  leads  to  the  conclusions  that 

\A/r=  .oHl  \J~Kj~ 

However,  this  treatment  presumes  that  the  conductivity 
a  is  temperature  independent.  In  fact,  it  varies  more 
than  an  order  of  magnitude,  e,g.,  for  copper,  the  con¬ 
ductivity  is  ~  10^  (ohm  cm)  ^  at  room  temperature  and 

~  it .  0  x  10^*  (ohm  cm)  at  ~  1500°K. 

Knoepfel  (p.  88)  obtains  an  approximate  solution 
to  the  "nonlinear"  (i.e.  ,  temperature  dependent  con¬ 
ductivity)  problem  for  the  case  where  the  magnetic 
field  is  strong.  He  writes  the  applied  field  as 

U(°)+)  -  ^ 


These  results  have  been  applied  to  the  example  of 
a  solid  current  carrying  copper  rod.  The  magnetic 
field  H  =  2I_=  1.256  (1/2-rrr)  =  1.256  (i/a).  In  Fig.  1, 
lOr 

the  loss  in  joules/cm  is  plotted  vs  time  for  a  variety 

of  values  of  I/a.  For  times  >  10  sec,  losses  can  be 
quite  large  in  a  practical  system  for  I/a  >.l  MA/cm. 

We  now  apply  these  results  to  the  design  of  a  sin¬ 
gle  turn  coil  of  radius  R  and  length  Z  inside  which  we 
wish  to  maximize  the  energy  density.  We  assume  the 
available  energy  is  1  megajoule  and  the  Inductance  L 

o 

of  the  feed  line  is  10  nh.  However,  for  purposes  of 
simplicity,  we  assume  also  that  the  width  of  the  feed 
line  is  so  great  that  losses  in  the  line  are  negligible, 
R  is  assumed  to  be  1  cm,  and  we  will  determine  the  val¬ 
ue  of  Z  which  maximizes  the  the  internal  energy  density 
of  the  coil. 

If  Z  »  R,  the  coil  inductance  is 

L  -  lb1  =  fle/ir/eS) 

c  /ooi  £ 


and  obtains  solution  valid  if  H  »  h  .  Later,  he  clar- 
ifies  the  interpretation:  hc  is  that  value  of  magnetic 

field  where  the  heating  due  to  the  energy  loss  has 
caused  the  conductivity  to  drop  by  a  factor  of  two. 
Thus,  if  the  applied  peak  magnetic  fields  are  generally 
below  h  ,  the  conductivity  variation  with  temperature 

is  probably  not  important. 

Since  our  interest  here  is  scaling,  we  adopt  the 
non-linear  treatment  of  Knoepfel.  The  revised  expres¬ 
sion  for  the  total  electromagnetic  energy  loss  is 

Wr  =  .  328  Hjy  ^  ) 

k 

The  calculations  leading  to  this  result  are  given  in 
Appendix  A. 

There  is  loss  mechanism  in  addition  to  joule  heat¬ 
ing  and  magnetic  field  penetration:  current  flowing 
in  the  conductor  generates  a  magnetic  field  outside  it; 
the  magnetic  pressure  pushes  on  the  conductor  surface 
and  generates  a  (weak)  shock  which  propagates  into  the 
conductor.  The  material  between  the  surface*the  shock 
front  is  compressed,  heated,  and  accelerated  by  the 
shock.  The  energy  required  to  do  this  must  come  from 
the  energy  source  and  thus  comprises  a  loss. 

The  treatment  of  this  problem  follows  Knoepfel, 
p.  257 >  and  the  calculations  leading  to  the  energy  loss 
per  unit  area  are  given  in  Appendix  B.  The  result  is, 
for  a  constant  applied  field, 


losses  are 

3 

* ;r=  • 

3z 

K 

H 

(Wt r%C 

c<  Aj 
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p  H1 


We  further  assume  H  is  constant  and  t  =  10  ^  sec. 

m 

The  energy  density  inside  the  coil  can  be  written  as 

B  =  -JI5I  Me*') 

V  TtPt-i 

If  Z  »  R,  H  =  .UttI/Z, 
and  I  can  be  obtained  from 

£0Z  l(LofLc)^  = 


The  expression  for  E/V  takes  the  form 

B  -  E<>  _  /k  j_  __  ga  j_ 

V  (LJZ+A,)*  ?3^  2rz 

Table  I  shows  how  E/V  varies  with  Z. 

TABLE  I 


Pq  is  the  normal  density  of  the  conductor,  and  S  and  C 
are  material  properties.  For  copper, 

S  =  1.5 

C  =  1+  x  10^  cm/sec  (~  speed  of  sound) 

P0  =  8.9  g/cm3 

For  very  large  currents,  the  magnetic  field  may  be 
large  enough  to  cause  the  second  term  in  the  square 
root  to  be  significant .  For  the  subsequent  analyses , 
however,  the  second  term  has  been  dropped,  and  in  this 
weak  field  approximation  the  shock  loss  per  unit  area 
is 

w  =  H  (Vs/ov^j 


Z(cm) 

1 

2 

3 

h 

5 

6 


E/V  (joules/cm^) 

-l.ll«10T 

-l+.l6xl05 

4.0x10^ 

1. 

6.7x10 

1* 

6.03x10 

A 

5.3x10 


Thus,  in  this  simplistic  calculation,  a  coil  length  of 
Z  ~  4.5  cm  gives  the  highest  internal  energy  density 

of  -  7x10^*  joules/cmk  corresponding  to  H  ~  2.9 
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Megagauss.  The  loss  is  10-15$,  and.  increases  rapidly  in  this  notation, 

if  R  or  L  is  decreased  or  if  E  is  increased, 
o  o 


Conclusions 


In  this  paper,  we  have  shown  that  conductors  car¬ 
rying  high  currents  show  two  kinds  of  losses:  electro- 

3 

magnetic  losses,  which  scale  as  H  ,  and  shock  com-- 

h  3 

pression  losses,  which  scale  as  H  or  H  ,  where  H  is 
the  self- field  at  the  conductor  surface.  Since  the 
loss  per  unit  of  conductor  area  depends  on  load  geome¬ 
try,  the  load  volume  cannot  he  made  arbitrarily  small 
to  achieve  the  highest  load  energy  density. 

In  a  vacuum  transmission  line  in  which  the  trans- 
r 

verse  electric  field  is  >  10  v/cm,  it  is  well  known 
that  a  current  carrying  sheath  forms  at  the  cathode 
and  about  half  the  current  flows  in  this  sheath.  If 
the  losses  in  the  sheath  are  negligible,  the  electro¬ 
magnetic  conductor  losses  In  the  cathode  conductor  in¬ 
ferred  from  the  discussion  of  this  paper  are  probably 
somewhat  high.  The  losses  at  the  anode  conductor  are 
unchanged.  Since  the  current  sheath  is  likely  to  be  in 
pressure  equilibirum  with  the  self  magnetic  field  and 
with  the  conductor  surface,  the  losses  due  to  shock 
compression  are  unaffected  by  the  presence  of  the 
sheath. 

Appendix  A:  Electromagnetic  Energy  Losses 

As  stated  in  the  test,  Poynting's  theorem  gives 


\A /_  -  -J-  jof 


tt(ot )  At 


W-r=  Sj.± 

T f  3 


Here,  I  is  an  integral  function  defined  by  Knoepfel  and 
and  T  (x)  is  the  gamma  function.  is  the  flux  pene¬ 

tration  skin  depth. 


fWe) 

P(X) 


(for  these  discussions,  Knoepfel 's  "n"  is  1.) 
Knoepfel  shows  that  the  ratios  W  /W^  and  are 


.kik 


.586 


Thus,  the  ohmic  and  magnetic  losses  are  nearly  equal. 
W„  is  given  by 


fir 


4. 

3 


PC Z) 


7 


.0^7  HzCojt) 


To  obtain  the  value  of  9H  x=o,  the  magnetic  diffusion 
3X 

equation  must  be  solved;  Knoepfel  (p.  53;  76-77)  ob¬ 
tains  for  constant  conductivity 


I 

2>  K  ~  Wo 


h0  pty)  ft  1  4  rr 

ro>  (rj  1 Ik>* 


where  H  (o,t)  =  H  fF~ 
0  V  T 


o 

The  total  energy  loss  per  unit  area 


is  com¬ 


posed  of  two  parts:  joule  heating  and  flux  penetration 
into  the  surface.  The  latter  is  a  loss  because  the 
flux  cannot  be  recovered  rapidly  from  the  skin  depth  of 
the  conductor.  Denoting  these  two  losses  by  and  W  , 
respectively, 


00 

WT  =  &rr  J  ^  (x,t^  dx 
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1 

00 

r  8H  (x ,t ’ ) 

1  2 

'U° 

11 

h> 

IS 

J  a, 

O 

J 

0 

l  3X 

dt 

These  results  apply  to  the  circumstance  where  0, 
the  material  conductivity,  is  constant.  However,  0  is 
actually  a  function  of  temperature.  For  high  currents 
and  large  magnetic  fields,  the  surface  layer  can  melt 
and  even  ablate,  and  the  conductivity  change  can  easily 
exceed  a  factor  of  two.  Knoepfel  (pp.  88-95)  treats 
the  "non-linear"  problem  (i.e.,  0  is  temperature-depen¬ 
dent)  for  the  case  where  the  surface  heating  is  a  domi¬ 
nant  effect.  He  starts  by  assuming 

tf(O)t)  -  y 

where  h  is  a  characteristic  field  for  the  material, 
c 

He  shows  that  h  is  that  magnetic  field  for  which  the 
c 

implied  surface  heating  causes  the  conductivity  to  drop 
by  just  a  factor  of  two.  For  purposes  of  scaling  to 
large  magnetic  fields  or  currents,  the  non-linear  re¬ 
sults  are  most  appropriate. 

Knoepfel 's  analysis  shows  how  to  modify  for  the 

non-linear  case.  He  estimates  in  two  ways  the  ratio  of 
the  non-linear  S  to  the  quantity 

<f> 

and  obtains 


Knoepfel  uses  his  solutions  to  the  diffusion  equation 
to  obtain  expressions  for  and  W  : 

*  hV°a)  5,  ■  *  rph)  ra)  i, 

wt  =  w\t)  .j_  P(3fx)  n (x)  z0 

Sir  3 


Kc 


!  «  ;  ^ 
J-J  Z  or  -i  l/i  -  .17!  Ctr  Jol 
3  ’  2. 
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We  will  use  the  average  of  the  numerical  factors, 
0-59-  Thus, 


WT  -  H  (°^)  4  s,  =  (ott)  4f 


ST 


m  H  ^  ^ 


So , 


y>  «.  * 

?TT 

for  which  the  solution  is 


a  -  £ 
<*s 


JTTW 

1  «n».c' 


This  value  will  he  used  for  further  calculations.  The 
division  of  this  energy  loss  between  joule  dissipation 
and  magnetic  losses  is  the  same  as  for  the  0  =  constant 
case. 

Appendix  B:  Shock  and  Compression  Losses 

The  treatment  of  this  problem  follows  Knoepfel, 
p.  257 5  and  starts  with  statements  of  conservation  of 
mass,  momentum,  and  energy: 


Pi  (Vs  -u)  =  PQ  Vs 
Pi  +  P,  (Vs  -u)? 


2  2 

The  quantity  H  S/(2irp  C  )  takes  on  the  value  1  if 
o 

H  ~  2-5  Megagauss.  When  H  «  2.5  Megagauss,  the  weak 
field  approximation  can  be  used: 


P  +  P  Vc 

o  o  G 


P  2  P 

E  +-i  +  -u  =  E  +~ 

1  Pi  2  o  Pq 

(see  Figure  2) 

Here,  p  =  density 

V-  =  shock  velocity 
b 

u  =  material  velocity 
P  =  pressure 

E  =  internal  energy  per  unit  mass 
Subscripts  0,1  refer  to  unshocked  and  shocked 
material,  respectively. 

We  assume  that  the  material  is  in  pressure  balance 
with  the  magnetic  field,  i.e. ,  (Fig.  2) 

P,  =  ^ 

1  sir 


From  the  momentum  conservation  equation  above, 
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when  H  >  2.5  Megagauss,  the  strong  field  approximation 
is  useful: 
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H_ 

\J  f-rrf.S 
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MS 

irptC1 
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The  work  done  by  the  piston  (the  magnetic  field)  is, 
for  a  constant  field. 


w. 


SH 


=  F,« t  ,Aur^a)t 
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For  the  strong  field  case. 


w<u  = 


rr 


H3t 
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In  the  weak  field  approximation, 


HgS 

2fp  C2 
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<<  1,  and 
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6«/  IfftC 


jS_  jr 
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Hence,  for  weak  fields  the  loss  scales  as  H**  (or  I/a)** 

3 

but  for  larger  fields,  the  scaling  goes  as  H  . 


8TT 


Since  P^  «  P  ,  and  from  the  mass  conservation  equa¬ 
tion, 

P,  -  *)  a  Vs  "  ^ 

O  ' 


Empirically,  it  is  known  that  a  linear  relationship  be¬ 
tween  Vt,  and  u  agrees  reasonably  well  with  reality; 

Vr,  =  C  +  Su,  C  and  S  are  material  constants 
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JOULES  /  cm 


ENERGY  LOSSES  IN  CIRCULAR  CONDUCTORS 


2 

Figure  1.  Losses  in  joules/cm  as  a  function  of  time  for  various  values  of  I/a  =  1/2ttx 
for  a  constant  current  in  a  round  solid  copper  rod.  The  solid  curve  is  the 
total  loss  (electromagnetic  plus  shock)  whereas  the  dashed  line  is  the 
electromagnetic  loss  only. 
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Figure  2.  Shock  compression  losses  in  a 

current  carrying  conductor.  The 
magnetic  field  is  produced  by  the 
current  flowing  in  the  conductor. 
The  quantities  p,  P,  and  u  refer 
to  the  material  density,  pressure 
and  velocity,  respectively;  the 
subscripts  o,l  refer  to  the 
unshocked  and  shocked  regions, 
respectively.  Vg  is  the  shock 
velocity. 
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